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Abstract. The phenomenological gravitational waveform models, which we refer to 
as PhenomA, PhenomB and PhenomC, generate full inspiral-merger-ringdown wave¬ 
forms of coalescing binary back holes (BBHs). These models are defined in the Fourier 
domain, thus can be used for fast matched filtering in the gravitational wave search. 
PhenomA has been developed for nonspinning BBH waveforms, while PhenomB and 
PhenomC were designed to model the waveforms of BBH systems with nonprecess- 
ing (aligned) spins, but can also be used for nonspinning systems. In this work, 
we study the validity of the phenomenological models for nonspinning BBH searches 
at low masses, mi ^2 > 4:Mq and mi -I- m 2 = M < 30Mq, with Advanced LIGO. 
As our complete signal waveform model, we adopt EOBNRv2 that is a time-domain 
inspiral-merger-ringdown waveform model. To investigate the search efficiency of the 
phenomenological template models, we calculate fitting factors by exploring overlap 
surfaces. We find that only PhenomC is valid to obtain the fitting factors better than 
0.97 in the mass range of M < 15M0. Above 15M0, PhenomA is most efficient in 
symmetric mass region, PhenomB is most efficient in highly asymmetric mass region, 
and PhenomC is most efficient in the intermediate region. Specifically, we propose an 
effective phenomenological template family that can be constructed by employing the 
phenomenological models in four subregions individually. We find that fitting factors 
of the effective templates are better than 0.97 in our entire mass region and mostly 
greater than 0.99. 

Keywords: gravitational waves, compact binary coalescence, phenomenological wave¬ 
forms 
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1. Introduction 

Coalescing binary black holes (BBHs) are among the most promising sources of 
gravitational wave (GW) transients for ground-based detectors, such as LIGO [1] and 
Virgo [ 2 ] . GW signals emitted from BBHs are conventionally divided into three phases: 
inspiral, merger and ringdown (IMR). Inspiral waveforms can be accurately obtained 
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by using post-Newtonian (PN) approximations. When the system reaches the ultra- 
relativistic regime, however, merger-ringdown waveforms should be calculated from 
numerical relativity (NR) simulations. A full IMR waveform can be constructed by 
combining a PN inspiral waveform and a NR merger-ringdown waveform. However, 
obtaining merger-ringdown waveforms from NR simulations is computationally very 
expensive. Thus, inspiral waveforms have mainly been used for low mass binaries below 
~ 30Mq in ground-based GW data analyses [3]. It has been anticipated that the 
inspiral waveforms would yield accurate analysis results for low mass binaries because 
in those systems the inspiral phase is likely to have most of the signal power in detection 
frequency band of ground-based detectors. 

However, several authors have pointed out that the absence of merger-ringdown 
phase can decrease the search efficiency, thus induce a non-negligible loss in detection 
rate even for low mass systems. When using inspiral-only template waveforms, a loss 
in signal-to-noise ratio (SNR) tends to increase with increasing total mass (M) of the 
signal, consequently the detection rate becomes reduced below 90% of that achievable 
by an IMR search if M exceeds some critical value Merit (this is explained in more 
detail at the end of section 3.1). Farr et al |1] obtained the values of M^it about 
10 — 15Mq using IMR signals generated from EOBNR model and TaylorF2 templates 
with the initial LIGO sensitivity. A similar study has been carried out by Brown et 
al [5] for Advanced LIGO, and they found M^it ~ IIAMq. Buonanno et al [6] 
also showed that M^it ~ 12Mq for various PN inspiral templates with EOBNR signals. 
Assuming phenomenological IMR signals, Ajith [7j and Gho [8] found that Merit ~ 15Mq 
for inspiral templates. While the former used 3.5PN TaylorTl templates, in the latter 
the inspiral waveforms were constructed by taking the inspiral parts into account from 
the original phenomenological IMR waveforms to avoid the systematic effect that arises 
when the model for signals is different from the model for templates. 

Efforts to establish the analytic IMR waveform models have been made over the 
past years. Those works have been carried out by constructing phenomenological 
families of waveforms in the Fourier domain by using PN-NR hybrid waveforms. The 
hrst phenomenological family (PhenomA), which can model the IMR waveforms of 
nonspinning BBHs, has been proposed by Ajith et al m El [ID]. The second family 
(PhenomB) has been developed by Ajith et al [H] by extending PhenomA to the case of 
nonprecessing spinning BBHs. The third family (PhenomG) has been made by Sanataria 
et al also for nonprecessing spinning BBHs. A recent phenomenological family 
(PhenomP) allows us to have IMR waveforms of processing BBHs [13]. On the other 
hand, time-domain IMR waveform models have also been developed using the effective- 
one-body (EOB) approach [TTl |T5]. Accuracies of EOB models have been improved 
by calibrating them to NR waveforms. The nonspinning EOBNR model described in 
Ref. [HI has been used in LIGO and Virgo to search for high mass BBHs for the hrst 
time ini and this model has been continuously upgraded as new NR simulations have 

f In the observational analysis of Ref. the sensitivity of a search with nonspinning EOBNR 
templates was estimated using simlulated PhenomA signals, which is roughly the opposite of our 
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been produced. The state-of-art calibrations of this model can be found in Refs. 

A nonprecessing spinning EOBNR model has also been proposed by Taracchini et al |20] 
and further calibrated IZH to 38 NR waveforms produced by the SXS Collaboration m 
Recently, Pan et al [23] have described a general procedure to generate precessing 
EOBNR waveforms starting from a nonprecessing spinning EOBNR model given in 
Ref. [20|, and they found remarkable agreement with two precessing NR waveforms. 
The EOBNR models of Refs. [IHl IMl El] have been implemented in the LSC Algorithm 
Library (LAL) [21] under the names of EOBNRv2, SEOBNRvl and SEOBNRv2. 

The stability of nonspinning EOBNR model has been investigated by Pan et al 
[22]. They showed that those waveforms of any length are sufficiently accurate for 
data analysis with advanced GW detectors. Recently, in addition, that has been 
anecdotally confirmed by Szilagyi et al [ 22 ] by using a long NR simulation in which the 
gravitational waveform was long enough to cover the entire frequency band of advanced 
GW detectors for a nonspinning binary at mass ratio mi/m 2 = 7 with a total mass as 
low as M = 45.5M0. They also found that existing phenomenological IMR waveforms 
display much greater disagreement with the NR simulation. EOBNRv2 is now believed 
to be sufficiently accurate to search for nonspinning BBHs with Advanced LIGO. 

However, EOBNRv2 waveforms are much slower to generate than the 
phenomenological models since a complicated system of ordinary differential equations 
has to be solved over a long time interval with small time steps. Moreover, time-domain 
waveforms need to be Fourier-transformed to perform the match, which is defined as 
a frequency-domain inner product weighted by the detector noise. We found that the 
speed of generating Fourier-transforming EOBNRv2 waveforms can be two orders of 
magnitude slower than the speed of generating phenomenological waveforms. Thus, the 
phenomenological waveform families can be used to densely cover the parameter space, 
particularly at low masses. However, the use of the phenomenological waveforms may 
result in a loss in search efficiency due to inaccuracies of the approximations. 

In this work, we adopt as our complete signal waveform model EOBNRv2 and test 
the validity of the phenomenological template models in searches for nonspinning BBHs 
with Advanced LIGO. We take into account all the spherical harmonics modes of the 
EOBNRv2 waveforms available in the LAL, those are the leading (2, 2) mode and the 
four sub-leading modes (2, 1), (3, 3), (4, 4) and (5 ,5). Because the phenomenological 
fitting coefficients of PhenomB are defined differently from those of PhenomA, a 
PhenomB waveform with masses (mi, m 2 ) and a zero spin is different from a PhenomA 
waveform with the same masses. A PhenomP waveform with masses (mi, m 2 ) and 
a zero spin is the same as a PhenomG waveform with the same masses and a zero 
spin because PhenomP is derived from PhenomG just by adding an effective precession 
spin. Therefore, we only consider PhenomA, PhenomB and PhenomG models for our 
nonspinning BBHs. When using PhenomB and PhenomG as our templates, we will 
neglect the spin effect by choosing zero spins in the wave functions. In this work, we 
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focus on low mass systems with mi, 2 > 4Mq and mi+m 2 = M < SOMq, and investigate 
the fitting factors for each template model. In particular, we show that for nonspinning 
BBH searches a template bank made up of several models is feasible and can be much 
more efficient. For this purpose, we propose an effective template family that is made 
up of the three phenomenological models and show that this template family can have 
the fitting factors better than 0.97 for all signals in our entire mass region. 


2. Phenomenological waveforms 

Obtaining phenomenological models involves finding fitting parameters of the analytic 
waveform family using PN-NR hybrid waveforms. The phenomenological waveform 
families are defined in the Fourier domain as the form, 

fcph=„(/)=.4,ff(/)e*'-»W, (1) 

where Aes{f) and Tefr(/) are the effective amplitude and the effective phase, and those 
are modeled separately. The model waveforms are parameterised by their physical 
parameters such as total mass M, symmetric mass ratio r] = 77111712 /M'^ and effective 
spin parameter y = (1 + (5)yi/2 + (1 — 5)x2/2 where 5 = (mi — m 2 )/M and Xi = Si/rrif, 
Si being the spin angular momentum of the ith BH. In this section, we only provide 
the functional forms of the phenomenological models. More details can be found in 
mum, [H] and |I2] for PhenomA, PhenomB and PhenomC, respectively. 


2.1. PhenomA 


Amplitude spectrum of a PhenomA waveform is divided by two phenomenological 
frequencies, /merg and /rmg, and terminates at the cutoff frequency fcut- 


(///merg)"^'^® if /< /merg 

(///merg)"^'^^ if /merg </< /ring 

W C{f, /ring, d) if /ring < / < /cut, 


Aeff = 71/-^/® 


' merg 


( 2 ) 


where A is the wave amplitude factor whose value depends on the binary masses and 
five extrinsic parameters determined by the sky location and the binary orientation. 
£(/, /ring, d) = if-f- '^) 2 +n 2/4 ^ Loreutzian function that has a width d, and that 

is centered around the frequency /ring- The normalization constant, w = ^ \T^) ’ 

is chosen so as to make Aefr(/) continuous across the transition frequency /ring- The 
parameter /merg is the frequency at which the power-law changes from /“^/® to 
The phenomenological parameters /merg, /ring, d and /cut can be obtained by fitting their 
formulas to the PN-NR hybrid waveforms, and those are finally given in terms of M 
and rj. The effective phase is expressed as 

'Feff(/) = 27r/tc + 0c +- ^{Xkrf+ yk'n + (3) 


where tc and 0c are the coalescence time and the coalescence phase. The coefficients 
introduced in the phenomenological parameters and the effective phase are tabulated 
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in Table 1 of [7]. The accuracy of PhenomA templates was examined by using (~ 30) 
PN-NR hybrid waveforms finely spaced in the parameter range mi/m 2 = O' < 4 and 
50 < M/Mq < 200, in which the fitting factors were found to be greater than 0.99 for 
Initial LIGO, Advanced LIGO and Virgo. 

2.2. PhenomB 

PhenomB corresponds to an extended version of PhenomA to nonprecessing spinning 
BBHs by incorporating the single spin parameter y. The effective amplitude is defined 
by 

All = .4/,’''“; Ml„, (///i)“^''“(l + if/l</</2 (4) 

[ Wr £(/, /2, d) if /2 < / < /cut, 

where v = (vrM/)^/^, 02,3 are given in terms of r] and x, ^ 1,2 are given in terms of x 
and the normalization constants Wm and Wr are chosen so as to make Aeff(/) continuous 
across the transition frequency /2 and /i, respectively. The effective phase is defined by 

4'efr(/) = 27r/G + 0e +(^1 +• (5) 

The coefficients introduced in this model are tabulated in table 1 of m- PhenomB 
templates were examined in the parameter range q < 10, and M < dOOM©, in which 
the fitting factors were greater than 0.965 for Initial LIGO. 


2.3. PhenomC 


PhenomG has also been developed for nonprecessing spinning BBHs. The wave 
amplitude terminates at /cut = 0.15/M, and that is constructed from two parts as 

Aeff = ApMif)wJ^ + Ard(/)w^, (6) 

where ApM is the premerger amplitude calculated by a PN inspiral amplitude with the 
addition of a higher order frequency term: 

Apm(/) = 44pn(/) + 7i/^^^, (7) 

V k=2 / 

where G = vrM/, and Ard is the ringdown amplitude: 

Ard = /RD(a, M), 52Q{a)]a)r'^/\ (9) 


where the Lorentzian function is defined by £'(/, /o, d) = d^/[(/ — foY + d/4], and Q is 
the quality factor which depends on the final BH spin a. The two amplitude parts can 
be combined by tanh-window functions: 



4(/ - /o) 


1 ± tanh 


d 


( 10 ) 
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Table 1. Parameter ranges of the phenomenological models valid for nonspinning 
BBH searches with the given detector sensitivity models. 


Model 

PhenomA 

PhenomB 

PhenomG 

Mass range [Mq] 

50 < M < 200 

M < 400 

M < 350 

Mass ratio range 

g < 4 

g < 10 

g < 4 

Detector 

Initial LIGO, Virgo, Advanced LIGO 

Initial LIGO 

Advanced LIGO 


where d = 0.005. The transition freqnency /o is determined by /o = 0.98 /r,d where /rd 
is a ringdown freqnency given in terms of M and a. The effective phase is calcnlated 
by a complete SPA inspiral phasing V'spa, a premerger phasing -^pm and a ringdown 
phasing as 

^efr(/) = '4’spawJ^ + V’PMta/jta/a + (11) 

with /i = 0.1/rd, /2 = /rd using d = 0.005 in the window functions. The premerger 
and ringdown phasing have the forms 

V'PM = ^ + a2f~^ + + 04 + 05 /^^^ + «6/), (12) 

h 

V’RD = /^l +/^2/, (13) 

where the oa: coefficients are inspired by the SPA phase, redefined and phenomenolog¬ 
ically htted to agree with the PN-NR hybrid waveforms, while /3i^2 parameters are not 
htted but obtained from the premerger ansatz by taking the value and slope of the phase 
at the transition point /rd- The coefficients introduced in this model are tabulated in 
table 2 of [ 12 ]. PhenomC templates were examined in the parameter range O' < 4, and 
M < 35 OM 0 , in which the model htted the original hybrid waveforms with the htting 
factors better than 0.97 for Advanced LIGO. 

2.4- Comparison 

In table we summarise the parameter ranges, in which the validity of the 
phenomenological models is conhrmed with the detector sensitivity models given in the 
last row. While the validity of PhenomA is tested only in high mass region, PhenomB 
and PhenomC are valid in broad mass region. PhenomB is applicable in wider range of 
mass ratio q < 10, while the others are valid only in the range of g < 4. 

Figure [T] shows the Fourier-domain amplitude spectra of EOBNRv2 and the 
phenomenological waveforms for a nonspinning binary with masses of (10,10)M©. 
The vertical lines indicate the transition frequencies that are calculated from 
phenomenological hts to the PN-NR hybrid waveforms. While PhenomA and PhenomB 
have two transition frequencies /merg and /ring, PhenomC has one frequency /q. In this 
hgure, we hnd that the phenomenological models show similar curves, and those agree 
well with the amplitude spectrum of the EOBNRv2 waveform. However, the phase 
rather than the amplitude is a main determining factor in the match calculation, so the 
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Frequency [Hz] 


Figure 1. Fourier-domain amplitude spectra of EOBNRv2 (E), PhenomA (A), 
PhenomB (B) and PhenomC (C) for a nonspinning BBH with masses (10,1O)M0. 
The transition frequencies /merg and /ring are denoted by the black (PhenomA) and 
red (PhenomB) lines, and /o (PhenomC) is denoted by the blue line. Small oscillations 
in the EOBNRv2 spectrum are due to an edge effect in the fast Pourier transform, and 
this does not affect our analysis. 


phase differences of the template waveforms from the EOBNR signals will mainly affect 
onr results of fitting factors. 


3. Validity of the phenomenological models and the effective template 
family 


3.1. BBH search and fitting factor 


The matched filter is the optimal filter for a signal of known shape in stationary 
Gaussian noise. Because the waveforms emitted from coalescing binaries can be modeled 
reasonably well, the matched filter can be employed in the search for BBHs. The match 
between a signal (hs) and a template (ht) is expressed using a standard inner product 
dehned by 


{hs\ht) = 4Re 


/low Snif) 


df, 


(14) 


where Sn{f) is a detector noise power spectral density and /low is a low frequency cutoff 
of the waveforms that depends on the shape of S'„(/). We use the zero-detuned, high- 
power noise power spectral density of Advanced LIGO izu, and assume /low = 10 Hz. 

When calculating a match, we use the normalized waveform h{f) = h{f)/{h\h)^^‘^, 
then the phenomenological models for nonspinning BBHs can be given in terms of 
M, 7], tc and fie for the match calculations. On the other hand, the match can be 













easily maximized over tc and 0c [2H], and the two mass parameters M and t] are main 
parameters in data analysis for nonspinning binary systems. Note that we use the chirp 
mass Me = Mrf^^ instead of M in this work. Therefore, the overlap is defined by 
maximizing the match over tc and 0 c as 


P 


max 


{hs\ht) 

\j {hs\hs)(ht\ht) 


(15) 


Finally, the fitting factor (FF) is the best-match between two waveforms maximized 
over all possible parameters [22] • In this work, thus, FF is obtained by maximizing the 
overlap over Me and rj as 


FF = maxP(Mc, r/). (16) 

Mc,r} 

For data analysis purposes, FF is used to evaluate the search efficiency. The gravitational 
wave searches use a bank of template waveforms constructed for the corresponding mass 
range [291 EDI EH 132]. Typically, a template bank requires that the total mismatch (i.e. 
1-FF) between the templates and signals does not exceed 3% [3lE3| including the effect 
of discreteness of the template spacing. In this work, in order to avoid the effect of 
discreteness, we use sufficiently fine spacings in the (Mc-t]) plan^ The detection rate 
is proportional to the cube of the SNR, and the SNR is calculated by p = (h 5 |hs)^/^FF. 
Thus, a FF = 0.97 corresponds to a loss of detection rates of ~ 10% . 


3.2. FFs of the phenomenological template models 


Since EOBNRv2 can model nonspinning BBH waveforms much more accurately 
compared to the phenomenological models, we employ EOBNRv2 as our complete signal 
model and the phenomenological models as template models. As described in Eq. (16), 
we calculate FF by exploring a two-dimensional overlap surface finely spaced in the (Mc- 
p) plane. In figure we summarise FFs of the phenomenological template models for 
nonspinning BBHs with masses mi, 2 > 4Mq and M < 3 OM 0 . The blue, red and black 
contours correspond to FF = 0.99, 0.98 and 0.97, respectively. This result shows the 
valid criteria of the phenomenological template models for nonspinning BBH searches 
at low mass region. Firstly, although PhenomA has been tested only in the high mass 
region M > 50Mq, our result shows that this model is also efficient in the region 
M > 12Mq. PhenomA has been tested in the region O' < 4, and FFs in our result are 
also sufficiently high above 0.97 in that region. Secondly, PhenomB has been examined 
for Initial LIGO in the range of M < 400Mq, but we find that FFs can be smaller 
than 0.97 for Advanced LIGO at the masses below ~ 15Mq. In addition, PhenomB has 
been examined in the range of q < 10, but our result shows that this model is invalid 
at symmetric mass region broadly. Finally, the PhenomG has been found to be valid in 


§ To obtain FF for one signal, for example, we repeat a grid search near the signal varying the search 
area and the template spacings until we can estimate a size of the overlap contour P — 0.995 roughly, 
where P is an overlap weighted by the maximum overlap value in that contour, and finally we find FF 
by performing a 51 x 51 grid search in the region P > 0.995. 
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m [M„] 


Figure 2. Fitting factors of the template models to EOBNRv2 signals for nonspinning 
BBHs. Darker shading indicates lower fitting factors. Dashed lines indicate constant 
mass ratios, mi/m 2 = 1,2,3,4 from top to bottom. For comparison, we also present 
the result of TaylorF2 inspiral template model. Note different mass ranges between 
TaylorF2 and the phenomenological models. 


broad mass region up to ~ SSOMq with q < 4 for Advanced LIGO, and our result also 
shows that FFs are overall better than 0.97 in the entire low mass region with q < 4. 

On the other hand, (Fourier-domain) inspiral-only templates can also be valid for 
BBH searches if masses of the signals are sufficiently small a El 0 m i. Brown et 
al [5] has investigated FFs of TaylorF2 inspiral templates to EOBNRv2 signals in the 
mass region 3Mq < mi,m 2 < 25Mq for Advanced LIGO (see, figure 3 therein). Their 
result showed that FFs tend to decrease with increasing total mass, and the trend of 
FFs is very simple compared to those of the phenomenological models. This is because 
TaylorF2 waveforms do not have MR phases and the cutoff frequencies are determined 
only by the total masses as nMfcut = 6“^/^. In figure we also present FFs of TaylorF2 
template model to EOBNRv2 signals only at the masses M < 16Mq. We find that the 
trend of FFs is very similar to that of Ref. |5] but our FFs are a bit better. A main 
reason of this difference is that Ref. |5] used less dense templates, consequently their 
result included the effect of discreteness of the template spacing. Our result show that 
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mi/m2=l 



M=30Mq 



ml m2 


Figure 3. Fitting factors calculated by changing a total mass with a fixed mass ratio 
mi/m 2 = 1 (left) and calculated by changing a mass ratio with a fixed total mass 
M = 30Mq (right). The dotted lines indicate M = IBMq (left) and q = 1.4,4 (right). 


TaylorF2 template model is valid (i.e. FF > 0.97) for nonspinning BBH searches if 
M < 15Mq, and FFs of that can be larger than 0.99 if M < 8Mq. However, a 
comparison of FFs between TaylorF2 and PhenomC shows that PhenomC is still more 
efficient than TaylorF2 even at very low masses in our mass region. 

3.3. Constructing the effective phenomenological template family 

In figure we can hnd that PhenomC is most efficient among the three models for 
nonspinning BBH searches although that has been developed for nonprecessing spinning 
BBHs. FFs of PhenomC are always better than 0.97 independently of the mass in the 
region g < 4. On the other hand, for an arbitrary signal model with masses (mi,m 2 ), 
we can choose one optimal model whose FF is larger than those of the other two models. 
The optimal htting factor {/F/F) is dehned by 

/F/F{mi,m 2 ) = Max[A(mi, m 2 ), B(mi, m 2 ), C(mi,m 2 )], (17) 

where a(mi,m 2 ) corresponds to a fitting factor calculated by using a template model 
a. Then, the optimal template family can be obtained by choosing /F/F for all signals 
in our mass region. 

More simply, if we divide the parameter region into four subregion and employ 
the models in the subregions individually, an effective template family can be easily 
obtained. The rule to divide the parameter region is motivated by the valid criteria of 
the models described in the previous subsection. In figure the left panel shows FFs 
calculated by changing a total mass with a hxed mass ratio g = 1, and the right panel 
shows FFs calculated by changing a mass ratio with a hxed total mass M = SOMq. In 
the left panel, only PhenomC is valid at the masses lower than ~ 15Mq, while PhenomA 
is most efficient at M > ISMq. In the right panel, PhenomA is most efficient at g < 1.4 
and PhenomB is most efficient at g > 4. In the region 1.4 < g < 4, all models have FFs 
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Figure 4. Fitting factors of the effective template family for nonspinning BBHs 
calculated by using PhenomA (A), PhenomB (B) and PhenomC (C) individually in 
the four subregions divided by the dashed lines. 


greater than 0.98, but the most efficient one is PhenomC. Note that g = 4 is consistent 
with the boundary of the valid criteria of PhenomA and PhenomC as described in 
table [Tj From these features, we divide the parameter region into four subregions and 
construct an effective phenomenological template family as in figure The effective 
template family is composed of PhenomC waveforms in the region M < 15Mq (as shown 
above PhenomC is overall better than TaylorF2 in this region, so TayorF2 is unnecessary 
for construction of the effective template family), while in the region M > 15Mq, that 
adopts individual models in three subregions divided by the mass ratios q = 1.4 and 4. 
As one can see, all FFs of the effective template family are better than 0.97 and mostly 
greater than 0.99. We found that the FF contours were almost consistent with those of 
the optimal template family. 

4. Summary and future work 

Next-generation ground-based GW detectors, such as Advanced LIGO [T] and advanced 
Virgo [2] , will start observations in coming years. Although astrophysical BHs are likely 
to have spins, nonspinning binaries may be the hrst search targets because nonspinning 
searches are computationally much cheaper than spinning searches. In addition, the 
Fourier-domain template models should also be used for fast matched filtering. In 
this context, we investigated the validity of the existing phenomenological template 
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models for nonspinning BBH searches at low masses with Advanced LIGO sensitivity 
assuming that EOBNRv2 can model the signal waveforms exactly. Although PhenomB 
and PhenomC have been developed for spinning systems, we also applied these models 
to nonspinning searches by choosing zero spins. We found that PhenomA is valid for 
low mass BBH searches in the region M > 12Mq and g < 4, and PhenomB is valid 
in the region M > 15Mq but invalid in broad asymmetric mass region. PhenomC is 
most efficient among the three models, with which FFs are always better than 0.97 in 
our mass region except highly asymmetric region. Specihcally, we proposed an effective 
phenomenological template family that could be easily constructed by employing the 
phenomenological models in the four subregions individually. The effective template 
family gives FFs better than 0.97 in our entire mass region and mostly greater than 
0.99. The effective template family can directly applied to the search analysis, thus that 
will allow us to conduct very efficient and fast analysis for the low mass, nonspinning 
BBH searches with forthcoming ground-based detectors. 

For spinning BBHs, PhenomB and PhenomC have been constructed by using 
nonprecessing PN-NR hybrid waveforms with several spin values mm- Recently, 
Piirrer [M] has also built a Fourier-domain reduced order model using time-domain 
nonprecessing spinning FOBNR waveforms (the state-of-art version of this model is 
SFOBNRv2ROM, which is implemented in the LAL). Using these three models, our 
study can be extended to nonprecessing spinning systems, and the feasibility of a 
three-dimensional effective template family can also be investigated when an exact 
hducial model for these systems is available. On the other hand, although we only 
considered the phenomenological models in this work, our approach on determining 
FFs can easily be applied to the SFOBNRv2ROM template model assuming zero spins, 
and our effective template family can be further improved by including the result of 
that model. We will compare our effective template family described in this work to the 
result of SFOBNRv2ROM, and investigate how much the effective template family can 
be improved by taking into account that model. 

The necessity of using full IMR waveforms is more pronounced in high mass region 
because a portion of merger-ringdown phase can signihcantly impact the total SNR 
(e.g., see hgure 9 in |31]). The phenomenological models are generally more efficient for 
high mass BBH searches. This will also be studied in detail in a future work. 
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